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1.  Introduction 

The  study  of  characteristic  functions  has  several  aspects.  Characteristic  func¬ 
tions  were  introduced  to  permit  the  application  of  powerful  analytical  methods 
in  probability  theory.  They  were  first  used  as  a  tool  to  study  limit  theorems, 
however  the  scope  of  their  applications  has  constantly  widened  and  includes 
now  a  large  variety  of  problems  in  probability  theory  and  in  mathematical 
statistics.  More  recently  mathematicians  began  to  investigate  problems  concern¬ 
ing  characteristic  functions  for  their  intrinsic  mathematical  interest.  In  some  of 
these  problems  their  probabilistic  or  statistical  origin  is  still  apparent,  in  others, 
the  analytical  character  becomes  dominant. 

The  present  survey  deals,  therefore,  with  a  variety  of  loosely  connected  topics. 
In  the  first  part,  we  study  certain  frequency  functions  whose  characteristic 
functions  are  known.  These  include  the  stable  distributions,  P61ya  type  distri¬ 
butions,  and  a  related  family.  These  results  are  interesting  from  the  probabilistic 
as  well  as  from  the  analytic  viewpoint.  The  second  part  deals  with  problems 
motivated  by  certain  statistical  questions.  The  third  part  treats  the  arithmetic 
of  distribution  functions  and  related  analytical  problems. 

Let  F{x)  be  a  distribution  function,  that  is  a  nonnegative,  right-continuous 
function  such  that  F(— <»)  =0  while  F(-|-<»)  =  1.  The  Fourier-Stieltjes  trans¬ 
form  of  F{x),  that  is  the  function 

(1.1)  /(<)  =  j^^e^dFix) 

is  called  the  characteristic  function  of  F{x).  In  this  paper  we  denote  distribution 
functions  by  capital  letters,  as  F{x),  and  the  characteristic  function  of  F{x)  by 
the  corresponding  small  letter,  asf{t).  If  subscripts  are  used  on  the  symbol  for  a 
distribution  function  then  the  same  subscripts  are  attached  to  its  characteristic 
function. 

This  work  was  supported  by  the  National  Science  Foundation  through  grants  NSF-G-4220 
and  NSF-G-9968. 
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Part  I.  Application  of  characteristic  functions  to  the 

STUDY  OF  DISTRIBUTION  FUNCTIONS 

It  is  often  of  interest  to  express  properties  of  distribution  functions  in  terms  of 
their  characteristic  functions.  In  this  part,  we  deal  mainly  with  stable  and  with 
unimodal  distributions.  However,  we  mention  here  three  interesting  isolated 
results.  J.  R.  Blum  and  M.  Rosenblatt  [1]  gave  conditions  which  assure  that  an 
infinitely  divisible  distribution  is  discrete  or  continuous  or  a  mixture.  J.  Shapiro 
[52]  derived  a  necessary  and  sufficient  condition  for  the  existence  of  moments 
of  an  infinitely  divisible  distribution.  E.  J.  G.  Pitman  [45]  studied  characteristic 
functions  which  have  a  derivative  of  odd  order  at  the  origin. 


2.  Stable  distributions  and  their  frequency  functions 

A  distribution  function  F(x)  is  said  to  be  stable  if  to  every  bi  >  0,  62  >  0,  Ci 
and  C2  there  corresponds  a  positive  number  h  and  a  real  c  such  that 

(2.1)  = 

holds  The  asterisk  denotes  here  the  operation  of  convolution.  Stable  distribu¬ 
tions  were  first  investigated  by  P.  L6vy,  who  showed  that  they  are  limit  laws  of 
normed  sums  of  independently  and  identically  distributed  random  variables. 
The  characteristic  function  of  a  stable  distribution  is  called  a  stable  characteristic 
function.  The  relation  (2.1)  can  be  expressed  in  terms  of  characteristic  functions 
as 

(2.2)  f(bit)f{b2t)  =  f(bt)e'^\ 
where  7  =  c  —  Ci  —  C2. 

It  is  possible  to  determine  all  stable  characteristic  functions  f{t);  they  are 
given  by 

(2.3)  log/a(<;  a,  |3,  c)  =  iat  -  c\t\^  |l  fa)(l<|,  a)|» 

where 


(2.4) 


a) 


.  Ta 

tan 


a  9^  1 


a  =  1 


and  where  c  ^  0,  [iS]  ^  1,  0  <  a  ^  2  while  a  is  a  real  number.  The  number  a 
is  called  the  exponent  of  the  stable  distribution. 

It  is  known  that  all  stable  distributions  are  absolutely  continuous;  we  denote 
the  frequency  functions  of  the  stable  distribution  with  parameters  a,  a,  j8,  c  by 
Pa(x]  a,  j8,  c)  and  write  p{x]  a,  /3,  c)  for  po{x;  a,  j8,  c)  and/(<;  a,  $,  c)  for  fo(t;  a,  c) . 
Then 
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(2.5) 


Vix)  a,  I3,c)  =  — 


a,  iS,  c)  dt. 


It  follows  easily  from  (2.3)  and  (2.5)  that 

(2.6)  p{x]  a,  c)  =  p(-x;  a,  -/3,  c). 

Explicit  expressions  for  the  frequency  functions  of  stable  distributions  are 
only  known  in  a  few  isolated  cases.  Thus  (2.5)  yields  for  o-  =  2  the  normal  fre¬ 
quency  function,  for  a  =  1  and  ^  =  0  the  density  function  of  the  Cauchy  dis¬ 
tribution.  P.  L4vy  and  N.  V.  Smirnov  determined  also  the  stable  distribution 
with  a  =  1/2;  it  belongs  to  the  system  of  Pearson  curves  (Type  V).  Apart  from 
these  three  cases,  no  stable  distributions  are  known  whose  frequency  functions 
are  elementary  functions. 

The  analytical  behavior  of  the  frequency  functions  of  stable  laws  was  first 
studied  by  A.  I.  Lapin  (see  [9],  p.  183),  who  showed  that  p(x]  a,  /3,  c)  is  an  entire 
function  if  a  >  1  while  for  a  =  1  the  radius  of  convergence  of  the  Taylor  series 
for  p{x;  a,  /3,  c)  in  the  neighborhood  of  a  point  x  of  the  real  axis  is  at  least  equal 
to  c.  The  case  a  <  1  was  investigated  by  A.  V.  Skorohod  [55]  who  obtained  the 
following  result. 

Theorem  2.1.  The  frequency  function  of  a  stable  distribution  with  characteristic 
exponent  a  <  1  has  the  form 


(2.7) 


p{x;  a,  c) 


~‘i>i(.T-“), 

X 


.T  >  0 

a:  <  0, 


where  ^>1(2)  and  <1*2(2)  are  entire  functions. 

V.  M.  Zolotarev  [68]  found  an  expression  for  the  density  of  a  stable  distribu¬ 
tion  with  exponent  a  greater  than  one  in  terms  of  a  density  with  exponent  1/a. 
In  view  of  (2.6)  it  is  sufficient  to  study  the  case  where  a;  >  0  while  —  1  ^  |8  ^  1 . 
Zolotarev  obtained  the  following  result. 

Theorem  2.2.  Let  v  =  tan  (7ra/2)  and  =  (1  +  then 


(2.8) 


avXp{—avX;  a,  /3,  1)  =  avX~^p 


) 


for  all  a;  >  0,  jjSl  ^  1  and  any  a  >  1.  Here  v  —  — j8  tan  (ira/2)  and 


(2.9) 


B  =  (an  tan  f-  arc  tan  v  -j-  (a 

L  2aJ  [^2a  L’r 


V.  M.  Zolotarev  [69]  studied  analytical  relations  between  stable  laws  with  dif¬ 
ferent  parameters  and  also  used  [70]  systematically  the  Mellin  transform  to 
investigate  analytical  properties  of  stable  laws.  In  the  first  of  these  two  papers 
he  considered  the  case  where  the  exponent  a  <  1  is  a  rational  number  and  showed 
that  p(x;  a,  /3,  1)  is  the  real  part  of  a  function  which  satisfies  an  ordinary  dif¬ 
ferential  equation.  Zolotarev  used  these  results,  and  some  known  formulas,  to 
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express  the  frequency  functions  of  stable  distributions  with  parameters  a  =  2/3, 
=  1;  a  =  3/2,  jS  =  1;  a  =  2/3,  /3  =  0;  a  =  1/3,  jS  =  1;  a  =  1/2,  jS  arbitrary, 
in  terms  of  higher  transcendental  functions. 

Somewhat  related  problems  were  treated  by  P.  Medgyessy  [43].  He  assumed 
that  the  exponent  a  —  m/n  is  a  rational  number  (w,  n  relatively  prime  integers) 
and  that  /3  =  0  in  case  a  =  1.  He  showed  that  p  =  p{x;  a,  /3,  c)  satisfies  a  linear 
partial  differential  equation  with  constant  coefficients. 


(2.10) 


=  0. 


Here  the  a,,  6i,  Ki,  with  i  =  1,  2,  3,  depend  on  m,  n  and  j8  but  are  not  uniquely 
determined.  This  fact  can  be  used  to  simplify  the  partial  differential  equation 
for  p.  In  case  a  9^  1  (but  not  necessarily  rational)  Medgyessy  [44]  obtained  two 
partial  integro-differential  equations  for  p.  Medgyessy  used  his  results  to  study 
the  decomposition  of  mixtures  of  stable  distributions. 

Among  the  an^ilytical  properties  of  stable  frequency  functions  that  were 
thoroughly  investigated,  we  must  finally  mention  their  asymptotic  behavior. 
A.  V.  Skorohod  [54]  gave  a  comprehensive  survey  of  asymptotic  formulas;  he 
considers  the  cases  shown  in  table  I.  More  recently  V.  M.  Zolotarev  [70]  gave 


TABLE  I 


a  <  1 

a  =  1 

a  >  1 

X  oc 

X  00 

X  — >•  «5 

/3  =  1 

X  — >  0(x  >  0) 

H 

i 

1 

8 

X  — >  — 00 

x—*x 

X  00 

j 

X  — ►  00 

-1  <  /3  <  1 

X  — 00 

X  — >■  — 00 

X  — »•  — 00 

X  — >  — 00 

X  — >•  — 00 

X  —00 

^  =  -1 

X  — »  0(x  <  0) 

X  — »  00 

X  — >  00 

asymptotic  expansions  for  stable  distributions  as  the  exponent  a  tends  to  the 
points  a  =  0  and  a  =  1. 

3.  Unimodality  of  certain  families  of  distribution  functions 

A  distribution  function  F{x)  is  said  to  be  unimodal  if  there  exists  at  least 
one  value  x  =  a  such  that  F{x)  is  convex  for  x  <  a  and  concave  for  x  >  a. 

A.  Wintner  [60],  [61]  was  one  of  the  first  to  investigate  unimodal  distributions. 
He  obtained 

Theorem  3.1.  The  limiting  distribution  of  a  sequence  of  unimodal  distributions 
is  unimodal. 
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Theorem  3.2.  T'he  convolution  of  two  symmetric  and  unimodal  distributions 
is  also  symmetric  and  unimodal. 

Another  important  result  is  due  to  A.  I.  Khinchin  [13]. 

Theorem  3.3.  The  function  f(t)  is  the  characteristic  function  of  a  unimodal 

distribution  if,  and  only  if,  it  can  he  represented  in  the  form  f{t)  =  {\/t)  g(u)  du 

where  g{u)  is  some  characteristic  function. 

The  study  of  unimodal  distributions  was  recently  stimulated  by  the  discovery 
of  an  error  in  a  proof  by  A.  I.  Lapin.  Lapin  asserted  that  the  convolutions  of 
unimodal  distributions  (not  only  of  symmetric  unimodal  distributions)  are 
unimodal.  Subsequently  Lapin’s  (erroneous)  statement  was  used  to  show  that 
the  distributions  of  the  L-class,  and  hence  also  all  stable  distributions,  are 
unimodal.  A  distribution  function  F(x),  and  its  characteristic  function /(O,  are 
said  to  belong  to  the  L-class  if  for  every  c,  with  0  <  c  <  1,  there  exists  a  charac¬ 
teristic  function  fc(t)  such  that  the  relation  f(t)  =  f(ct)fe(t)  holds  for  all  t.  The 
distributions  (characteristic  functions)  of  the  L-class  also  are  called  self-decom¬ 
posable  distributions  (characteristic  functions).  The  error  in  Lapin’s  proof 
was  pointed  out  by  K.  L.  Chung  [3]  (see  also  appendix  II  of  Gnedenko  and 
Kolmogorov’s  book  [9])  who  gave  also  counterexamples.  The  question  whether 
the  self-decomposable  and  stable  distributions  are  unimodal  became  therefore 
once  more  an  open  problem. 

1.  A.  Ibragimov  [10]  introduced  the  concept  of  strongly  unimodal  distribu¬ 
tions.  He  calls  a  distribution  strongly  unimodal  if  its  convolution  with  any  uni¬ 
modal  distribution  is  unimodal.  Ibragimov  showed  that  a  (nondegenerate)  dis¬ 
tribution  function  F(x)  is  strongly  unimodal  if,  and  only  if,  F{x)  is  continuous 
and  the  function  log  F'{x)  is  concave  on  the  set  of  points  on  which  neither  the 
right  nor  the  left  derivative  of  F(x)  vanishes. 

I.  A,  Ibragimov  [11]  constructed  an  example  of  a  characteristic  function  of 
the  L-class  which  does  not  belong  to  a  unimodal  distribution. 

A.  Wintner  [60]  showed  that  all  symmetric  stable  distributions  are  unimodal, 
he  [62]  also  succeeded  in  proving  the  unimodality  of  the  symmetric  and  self- 
decomposable  distributions.  Finally  I.  A.  Ibragimov  and  K.  E.  Cernin  [12] 
showed  that  all  stable  distributions  are  unimodal.  It  would  be  interesting  to 
obtain  a  characterization  of  all  unimodal  distributions  of  the  L-class. 

R.  G.  Laha  [18]  recently  obtained  two  sufficient  conditions  which  insure  that 
a  real  valued,  even  function  f(t)  is  the  characteristic  function  of  a  unimodal 
distribution.  These  conditions  are  similar  to  P61ya’s  condition  [46]  and  are  easily 
applicable. 

Theorem  3.4.  Let  f(t)  he  a  real  valued  and  even  function  which  satisfies  the 
conditions 

(i)  m  =  1, 

(ii)  limifi^»/(0  =  0, 

(iii)  the  function  tf'{t)  exists  and  is  continuous  for  all  real  t,  moreover  limt_+o(f ’(0 
=  limiti->oo^/'(0  =  0, 
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(iv)  the  function  g{t)  =  f{t)  +  tf{t)  is  convex  for  <  >  0. 

Then  fit)  is  the  characteristic  function  of  a  symmetric  and  unimodal  distribution. 

The  function  g{t)  satisfies  the  conditions  of  P61ya’s  theorem  (see  [46],  theorem 
1)  and  is  therefore  a  characteristic  function.  It  follows  from  the  definition  of  g{t) 

that  fit)  =  (1/0  giu)  du;  according  to  theorem  3.3  the  function  fit)  is  then 

the  characteristic  function  of  a  unimodal  distribution. 

Theorem  3.5.  Let  fit)  he  a  real  valued,  continuous  and  even  function  of  the 
real  variable  t  such  that  fiO)  =  1.  Suppose  that  there  exists  a  function  Aiz)  of  the 
complex  variable  z,  where  z  =  t  iy  with  t,  y  real,  that  satisfies  the  conditions 

(i)  /(O  =  ^  it)  for  real  t  >  0, 

(ii)  Aiz)  is  regular  in  the  region  —  ei  <  arg  2  <  7r/2  +  62,  61  >  0,  62  >  0, 

(iii)  1A(2)|  =  0(1)  as  \z\  -^0,  1A(2)1  =  Oi\z\~“)  as  I2I  — ^  00  with  a  >  1, 

(iv)  ImAiiy)  ^  0  for  y  >  0, 

where  Aiz)  may  also  assume  complex  values.  Then  fit)  is  the  characteristic  function 
of  an  absolutely  continuous,  symmetric  and  unimodal  distribution. 

It  follows  from  assumption  (iii)  that/(0  is  absolutely  integrable;  the  integral 

Q-iix  f(j^  di  =  -  /  cos  ixfit)  dl 

Jo 

exists  and  is  a  real  valued,  continuous  and  even  function  of  x.  We  use  assumption 
(i)  to  write  (3.1)  in  the  form 


(3.1) 


(3.2) 


pix)  =  -  Re  /"  e^^^  fit)  df  =  -  Re 
T  jo 


/: 


e^^^  Ait)  dt. 


To  evaluate  the  last  integral  we  consider  the  integral 


(3.3) 


along  a  contour  C  which  consists  of  the  following  four  parts:  (1)  the  segment 
r  ^  t  ^  R  of  the  real  axis,  (2)  the  circular  arc  T  of  radius  R  with  center  at  the 
origin  located  in  the  first  quadrant,  (3)  the  segment  r  Sy  ^  R  oi  the  imaginary 
axis,  (4)  the  circular  arc  7  of  radius  r  with  center  at  the  origin  located  in  the 
first  quadrant.  We  see  from  assumption  (ii)  that 

(3.4)  j  6^=^ Aiz)  dz  =  0, 

c 

we  also  deduce  from  (iii)  that,  for  a:  ^  0, 

(3.5)  lim  f  e’**A(z)  dz  =  lim  f  e^’^Aiz)  dz  =  0, 

T — ^0  »  R — ^  * 

T  r 


so  that 
(3.6) 


U  e‘‘'A(()  dt  =  i  e~^A(iy)  dy. 
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We  conclude  therefore  from  (3.2),  (3.6)  and  assumption  (iv)  that 

(3.7)  p{x)  =  —-  f  e-^^lm[A{iy)]dy, 

’’’  Jo 

so  that  p{x)  ^  0  for  all  real  x.  Let  0  <  rci  <  x^,  it  follows  from  (3.7)  and  (iv)  that 

(3.8)  p(0)  >  p(xi)  >  p{x^, 

hence  p{x)  has  a  unique  maximum  at  x  =0.  This  completes  the  proof  of  the 
theorem. 

The  application  of  these  theorems  yields 
Theorem  3.6.  The  function 

(3.9)  m  = 

is,  for  any  a  in  the  interval  0  <  a  S  2,  the  characteristic  function  of  a  unimodal 
distribution . 

The  fact  that  the  function  f{t),  as  dehned  in  theorem  3.6  is  a  characteristic 
function  was  already  established  by  Yu.  V.  Linnik  [26],  If  0  <  a  ^  1,  theorem 
3.6  follows  from  theorem  3.4;  if  1  <  a  <  2  we  obtain  it  from  theorem  3.5.  For 
a  =  2  formula  (3.9)  yields  the  characteristic  function  of  the  Laplace  distribution 
which  is  known  to  be  unimodal.  It  is  easy  to  show  that  f(t),  as  given  by  (3.9), 
cannot  be  a  characteristic  function  ifa<0orifQ:>2. 

A.  Wintner’s  [60]  result,  that  all  symmetric  stable  distributions  are  unimodal, 
follows  easily  from  theorem  3.6.  Let  w  be  a  positive  integer  and  write  On  = 

The  function 


(3.10)  gr.it)  =  [/(a„0]-"  = 

is,  according  to  theorem  3.2,  the  characteristic  function  of  a  symmetric,  uni¬ 
modal  distribution.  We  see  from  theorem  3.1  that  the  same  statement  holds  for 
the  function 

(3.11)  lim  g„it)  = 


Part  II.  Characterization'  problems 

Let  Xi,  X2,  •  •  •  ,  X„  be  a  sample  from  a  given  population  or,  more  generally, 
n  independently  but  not  necessarily  identically  distributed  random  variables. 
For  the  sake  of  convenience  we  shall  use  in  the  following  the  term  statistic  for 
a  single  valued  and  measurable  function  of  the  Xi,  X 2,  •  •  •  ,  Yn  even  in  case 
these  random  variables  are  not  identically  distributed.  In  this  part  we  consider 
the  problem  of  characterizing  the  distributions  of  the  random  variables  Xi, 
X2,  •  •  •  ,  X„  by  properties  of  certain  statistics.  The  work  done  in  this  direction 
up  to  1955  was  surveyed  at  the  Third  Berkeley  Symposium  [38].  However, 
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research  on  characterization  problems  continued  during  the  past  five  years.  As 
a  result  of  these  recent  investigations  the  problem  has  changed  its  character. 
Up  to  1950  one  studied  in  this  connection  mainly  a  few  specific  populations, 
while  the  more  recent  developments  deal  primarily  with  the  analytical  properties 
of  the  characteristic  functions  of  the  random  variables  Xi,  Xz,  •  •  •  ,  X„. 

In  section  4,  we  report  on  some  developments  which  are  still  closely  connected 
with  the  earlier  work.  In  section  5,  we  treat  some  of  the  recent  results  which  have 
a  more  analytical  character. 

4.  Characterization  of  certain  populations 

In  this  section  we  characterize  populations  by  means  of  certain  regression 
properties.  Let  Y  and  X  be  two  random  variables  and  assume  that  the  first 
moment  of  Y  and  the  kth  moment  of  X  exist  and  write  E(Y\X)  for  the  condi¬ 
tional  expectation  of  Y,  given  X.  We  introduce 

Definition  4.1.  The  random  variable  Y  is  said  to  have  pohjnomial  regression 
of  order  k  on  X  if  the  relation 

(4.1)  E{Y\X)  =  ^0  +  iSiX  +  •  • .  +  ^,X^  ^  0, 

holds  almost  everywhere. 

If  =  2  we  use  the  term  quadratic  regression,  if  it  =  1  we  speak  about  linear 
regression,  li  k  =  0,  that  is  if  E{Y\X)  =  E{Y)  almost  everywhere,  then  we  say 
that  F  has  constant  regression  on  X.  The  following  lemma  is  a  simple  generaliza¬ 
tion  of  a  result  (lemma  6.1)  in  [38]  and  is  proven  in  the  same  manner. 

Lemma  4.1.  Let  X  and  Y  he  two  random  variables  and  assume  that  the  expecta¬ 
tions  E(Y)  and  E(X’‘)  exist  where  k  is  a  nonnegative  integer.  The  random  variable 
Y  has  polynomial  regression  of  order  k  on  X  if,  and  only  if,  the  relation 

(4.2)  ECYe^*^)  =  £  /S,-  E{X’e^^^) 

j  =  0 

holds  for  all  real  t. 

In  an  earlier  paper  [38]  the  author  characterized  the  Poisson  population  by 
making  the  following  two  assumptions:  (i)  the  population  distribution  function 
has  the  point  x  =  0  as  its  left  extremity,  (ii)  the  statistic  S  =  kp+z  —  kp  has 
constant  regression  on  ki.  Here  p  ^  1  is  a  positive  integer  and  kj  denotes  the 
fc-statistic  of  order  j,  a  symmetric  and  homogeneous  polynomial  statistic  of  order 
j  such  that  E(kj)  =  k,  where  kj  is  thejth  cumulant  of  the  population  distribution 
function.  This  result  was  recently  generalized  [40]  by  considering  instead  of  S 
the  difference  of  two  ^-statistics  of  arbitrary  order  and  by  removing  the  restric¬ 
tion  that  the  population  was  one  sided.  The  following  result  was  obtained. 

Theorem  4.1.  Let  Xi,  Xz,  •  •  •  ,  Xn  be  a  sample  from  a  population  with  popu¬ 
lation  distribution  function  F{x)  and  denote  by  p  ^  1,  r  ^  1  two  positive  integers. 
Assume  that  the  moment  of  order  p  r  of  F(x)  exists.  The  population  distribution 
F{x)  has  the  characteristic  functions 
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(4.3)  f(t)  =  exp  -  1)  +  6rX2(e“*'  —  1)  +  Ciit  +  | 

if,  and  only  if,  the  statistic  kp+r  —  kp  has  constant  regression  on  k\.  Here 
Xi,  X2,  Cl,  C2  are  real  constants  such  that  Xi  ^  0,  X2  ^  0,  C2  ^  0  while  8r  =  [?'/2]  — 
[(r  —  l)/2]  and  €p  =  1  if  p  >  1  and  ei  —  0. 

The  symbol  [a;]  denotes  the  largest  integer  not  exceeding  x. 

This  theorem  characterizes  a  family  of  distributions  which  consists  of  con¬ 
volutions  of  a  Poisson  distribution,  the  conjugate  to  a  Poisson  distribution  and 
a  normal  distribution.  For  certain  values  of  r  and  p  one  or  two  of  these  factors 
must  be  absent.  This  theorem  shows  that  the  assumption  that  kp+r  —  kp  has 
constant  regression  on  ki  is  not  sufficient  to  characterize  the  Poisson  population. 
If  one  wishes  to  characterize  the  Poisson  population  one  must  impose  some 
additional  restriction. 

Theorem  4.2.  Let  Xi,  X^,  •  •  •  ,  X„  be  a  sample  from  a  population  with  popu¬ 
lation  distribution  function  F(x)  and  denote  by  p  1,  r  1  two  positive  integers. 
Assume  that 

(i)  the  {p  -f-  r)th  moment  of  F{x)  exists 

(ii)  F{x)  =  0  for  a:  <  0  while  F(x)  >  0  for  a:  ^  0. 

The  population  is  a  Poisson  population  if,  and  only  if  kp+r  —  kp  has  constant 
regression  on  k\. 

For  the  proof  of  these  results  certain  extensions  of  the  theorem  of  Marcinkiewicz 
are  needed  (see  [39]).  These  provide  necessary  conditions  which  an  entire  func¬ 
tion  must  satisfy  in  order  to  be  a  characteristic  function. 

M.  C.  K.  Tweedie  [59]  investigated  the  regression  of  the  sample  variance  on 
the  sample  mean.  He  characterized  several  populations  by  assuming  that  the 
sample  variance  has  quadratic  or  linear  regression  on  the  sample  mean.  More 
recently,  E.  G.  Laha  and  E.  Lukacs  [19]  considered  a  general  quadratic  statistic 

(4.4)  Q  =  Z  Z  aijXiXj  -h  i:  bjXj 

i-i }=i  y=i 

and  determined  all  populations  which  have  the  property  that  Q  has  quadratic 
regression  on  the  statistic  Xi  -f-  X2  -f  •  •  •  +  X„.  In  this  study  it  is  necessary  to 
distinguish  several  cases  which  are  defined  in  terms  of  relations  between  the 
coefficients  a^  and  by  of  Q  and  the  regression  coefficients  /Jo,  i8i,  The  following 
distributions  were  obtained:  (i)  the  normal  distribution,  (ii)  Poisson-type  dis¬ 
tributions  (that  is,  Poisson  distributions  with  scale  and  location  parameters), 

(iii)  binomial  and  negative-binomial  distributions,  (iv)  Gamma  distributions, 
(v)  distributions  with  characteristic  function 

(4.5)  f{t)  =  e*''‘[cosh  at  -b  i\  sinh  at\~% 

where  a,  X,  y,  p  are  real  and  a  5^  0,  p  >  0. 

We  conclude  this  section  by  mentioning  two  characterizations  of  the  Wiener 
process.  V.  P.  Skitovich  [53]  obtained  the  following  result. 

Theorem  4.3.  Let  X(t)  be  a  homogeneous  stochastic  process  with  independent 
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increments  which  is  defined  in  a  closed  interval  {_A,  B].  Let  a{t)  and  h(t)  he  two 
functions  which  are  continous  in  [A,  B]  such  that  I  [a{t)h{t)]^  dt  9^  0  and  that  at 
least  one  of  the  integrals 


(4.6) 


/: 


aKt) 

hKt) 


dt  or 


dt 


exists.  Suppose  that  the  two  stochastic  integrals  a{t)  dX{t)  and  b{t)  dX(t) 

are  independently  distributed,  then  X{t)  is  a  Wiener  process. 

The  stochastic  integrals  are  here  stochastic  limits  of  Riemann-Stieltjes  sums; 
however,  the  mode  of  stochastic  convergence  differs  somewhat  from  the  usual 
definitions.  V.  P.  Skitovich  requires  that  the  distribution  function  of  the  Riemann 
sums  should  converge  uniformly  to  the  distribution  of  a  random  variable.  This 
random  variable  is  then  called  the  stochastic  integral.  V.  P.  Skitovich  also  ex¬ 
tended  his  earlier  result  concerning  the  independence  of  linear  forms  to  linear 
forms  in  infinitely  many,  identically  distributed  random  variables. 

R.  G.  Laha  and  E.  Lukacs  [20]  characterized  the  Wiener  process  by  means  of 
the  following  regression  property. 

Theorem  4.4.  Let  X{t)  be  a  stochastic  process  which  satisfies  the  conditions 

(a)  X{f)  is  defined  in  a  finite,  closed  interval  [A,  B], 

(b)  X{t)  is  homogeneous  and  has  independent  increments, 

(c)  X(t)  is  of  second  order  and  its  mean  value  function  and  covariance  function 
are  of  bounded  variation  in  [A,  Bf 

Suppose  that  a{t)  and  h{t)  are  two  continuous  functions  defined  in  [A,  B]  such 
that  ait)h(f)  9^  0  for  all  t  E  [Ai,  Bi],  where  A  ^  Ai  <  Bi  ^  B.  Suppose  further 
that  a{t)  is  not  proportional  to  b{t).  Let 

(4.7)  Y  =  a{t)  dX{t),  and  Z  =  b{t)  dX{t) 

be  two  stochastic  integrals,  defined  as  limits  in  the  mean.  The  process  X(t)  is  a 
Wiener  process  if,  and  only  if, 

(i)  Y  has  linear  regression  on  Z 

(ii)  The  conditional  variance  of  Y,  given  Z,  does  not  depend  on  Z. 


5.  Analytical  aspects  of  the  characterization  problem 

The  characterization  problems  which  we  studied  in  the  preceding  section,  and 
also  those  treated  in  part  II  of  [38],  dealt  with  specific  statistics,  such  as  the 
sample  mean  and  A:-statistics.  Projjerties  of  the  distributions  of  these  statistics 
were  used  to  characterize  completely  certain  populations,  except  for  the  numeri¬ 
cal  values  of  some  parameters.  The  solution  of  these  problems  usually  was  carried 
out  in  three  steps:  (1)  The  assumptions  concerning  the  distribution  of  the  statis¬ 
tics  are  used  to  derive  a  differential  equation  for  the  characteristic  function  of 
the  population  or  the  characteristic  functions  of  the  random  variables  in  case 
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they  are  not  identically  distributed,  (2)  this  differential  equation  is  solved,  (3) 
those  solutions  are  determined  that  are  characteristic  functions.  Frequently  the 
last  step  is  the  most  difficult,  sometimes  the  differential  equation  is  so  compli¬ 
cated  that  it  cannot  be  solved  readily.  These  difficulties  lead  then  to  investiga¬ 
tions  concerning  the  analytical  properties  of  the  solutions  of  those  differential 
equations  that  are  characteristic  functions. 

The  first  results  in  this  direction  were  obtained  by  A.  A.  Zinger  and  Yu.  V. 
Linnik  [66],  [67].  To  formulate  these  results  it  is  necessary  to  introduce  the 
following  terminology. 

We  consider  an  ordinary  differential  equation 

(5.1)  Z  •  •  •/(^“>«)  =  c[/(0]" 

and  denote  the  order  of  the  differential  equation  by  m.  The  are  real  con¬ 

stants  and  the  sum  is  taken  over  all  nonnegative  integers  which  satisfy  the  con¬ 
dition 

(5.2)  il  +  i2  +  •  •  •  +  jn  ^  p. 

Here  p  is  an  integer  such  that  at  least  one  of  the  coefficients  with  ji  + 

•  ’  •  +  in  =  p  is  different  from  zero.  We  adjoin  to  the  differential  equation  (5.1) 
a  polynomial 

(5.3)  A(xi,  a:2,  •  •  •  ,  rtn)  =  ^  Z*  Z  Aj,..  •  ‘xt, 

where  the  first  summation  Z*  runs  over  all  permutations  {ki,  A:2,  •  •  •  ,  kn)  of  the 
first  n  positive  integers  while  the  second  summation  is  taken  over  all  subscripts 
satisfying  (5.2). 

The  differential  equation  (5.1)  is  said  to  be  positive  definite  if  the  adjoint 
polynomial  (5.3)  is  nonnegative. 

A.  A.  Zinger  and  Yu.  V.  Linnik  obtained 

Theorem  5.1a.  Suppose  that  the  function  f{t)  is,  in  a  certain  neighborhood  of 
the  origin,  a  solution  of  the  positive  definite  differential  equation  (5.1)  and  assume 
that  m  n  —  1.  If  the  solution  is  a  characteristic  function ,  then  it  is  necessarily 
an  entire  function. 

We  indicate  here  also  the  motivation  for  this  study.  Let 

(5.4)  P  =  P(Xi,  X2,  •  •  •  ,  X„)  =  Z  Aj,. .  .jXl  ■  •  Ai- 

be  a  polynomial  statistic  of  degree  p.  We  adjoin  to  P  the  statistic 

(5.5)  f  L*  L  ■  -n, 

where  the  summations  Z*  and  Z  are  taken  in  the  same  way  as  in  formula  (5.3). 

The  statistic  P  is  said  to  be  a  regular  polynomial  statistic  of  degree  p  and  order 
m  if  the  following  three  conditions  are  satisfied: 

(i)  The  statistic  P*,  adjoint  to  P,  is  a  nonnegative  polynomial, 

(ii)  no  exponent  in  P  exceeds  m, 
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(iii)  at  least  one  variable  in  P  has  the  exponent  m. 

The  assumption  that  a  regular  polynomial  statistic  P  has  constant  regression 
on  the  sample  mean  X  =  (Xi  +  X2  +  •  •  •  4-  Xn)/n  leads  then  to  a  positive 
definite  differential  equation  for  the  characteristic  function  of  the  population. 
In  this  case,  one  can  apply  theorem  5.1a  and  obtain 

Theorem  5.1b.  Let  Xi,  X^,  •  •  •  ,  Xn  be  a  sample  from  a  population  with  dis¬ 
tribution  function  F(x)  and  assume  that  F{x)  has  moments  up  to  order  m.  Let  A  be 
the  sum  Xi  +  X2  +  •  •  •  +  X„  and  let 

(5.6)  P  = 

be  a  regular  polynomial  statistic  of  degree  p  and  order  m,  where  m  ^  p.  If 

(i)  P  has  constant  regression  on  A, 

(ii)  m  ^  n  —  I, 

then  the  characteristic  function  f(t)  of  F{x)  is  an  entire  function. 

Theorem  5.1a  is  a  very  interesting  result  concerning  the  analytic  properties 
of  the  solutions  of  certain  ordinary  differential  equations.  If  one  wishes  to  use 
it  in  connection  with  characterization  problems  one  obtains  theorem  5.1b. 
However,  in  applying  theorem  5.1b  one  is  greatly  handicapped  by  the  severe 
restrictions  contained  in  its  assumptions.  This  is  illustrated  by  the  following 
facts.  The  positive  definiteness  of  the  differential  equation  excludes  the  /b-sta- 
tistics  kp  of  order  p  >  2,  whereas  we  know  that  the  normal  distribution  is 
characterized  by  the  property  that  kp  has  constant  regression  on  A  for  any  p  ^  2. 
Moreover,  in  the  case  of  ^2,  condition  (ii)  of  theorem  5.1b  restricts  the  sample 
size  n  to  3.  It  would  therefore  be  desirable  to  derive  similar  results  under  modified 
conditions. 

In  the  early  attempts  to  characterize  populations  by  the  independence  of  two 
statistics  it  was  always  assumed  that  certain  moments  of  the  population  dis¬ 
tribution  function  exist.  Later  it  was  possible  to  relax  or  remove  this  assumption 
in  special  cases.  For  example,  K.  C.  Chanda  [2]  and  Yu.  V.  Linnik  [29]  study 
the  independence  of  a  polynomial  and  a  linear  statistic  and  show  that  it  is 
sufficient  to  suppose  that  the  moment  of  order  5  >  0  exists,  where  5  is  not 
necessarily  an  integer  and  can  be  arbitrarily  small. 

A.  A.  Zinger  [63],  [64]  investigated  the  independence  of  two  polynomial 
statistics  and  eliminated  the  assumption  concerning  the  existence  of  moments. 
For  the  formulation  of  his  results  we  need  the  following  definition. 

A  polynomial  P  =  P(xi,  X2,  •  "  ,  Xn)  of  degree  m  is  said  to  be  admissible,  if 
the  coefficients  of  the  terms  xf,  where  j  =  1,  2,  •  •  >  ,n,  are  not  zero.  Here  and  in 
the  following  we  assume  that  similar  terms  have  been  collected  in  every  poly¬ 
nomial  which  we  consider. 

Theorem  5.2.  Let  Xi,  X2,  •  •  *  ,  Xn  be  n  independently  (but  not  necessarily 
identically)  distributed  random  variables.  Let  Pi{Xi,  •  •  •  ,  Xn)  and  P2(Xi,  •  •  •  ,  X„) 
be  two  admissible  polynomials.  If  Pi  and  P2  are  independent,  then  each  Xj,  where 
j  =  1,  2,  •  •  •  ,n,  has  finite  moments  of  all  orders. 

A  stronger  result  is  obtained  if  one  of  the  polynomials  is  a  linear  statistic. 
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Theorem  5.3.  Let  X\,  X2,  •  •  •  ,  X«  be  n  independently  {but  not  necessarily 
identically)  distributed  random  variables  with  charcuiteristic  functions  j\{t),f2{t),  •  •  •  , 
fn{t)  respectively.  Let  P  =  P{Xi,  X2,  •  •  •  ,  XJ  be  an  admissible  polynomial  sta¬ 
tistic  and  A  =  A(Xi,  X2,  •  •  •  ,  XJ  =  ajXj,  with  aj  ^  0  forj  =  1,  2,  •  •  •  ,  n, 
be  a  linear  form.  If  P  and  A  are  independent  then  the  characteristic  functions  fj{t), 
for  j  =  1,  2,  •  •  •  ,  n,  are  entire  functions  of  finite  order. 

Theorems  5.2  and  5.3  are  due  to  A.  A.  Zinger  who  showed  that  they  are  valid 
also  for  admissible  quasipolynomial  statistics.  A  function  S  =  22,  •  •  •  ,  Xn) 

is  called  a  quasipolynomial  if  there  exists  a  continuous  function  T{x)  and  two 
nonnegative  polynomials  Pi  =  Pi{xi,  0:2,  •  •  •  ,  Xn)  and  P2  =  P2{xi,  X2,  •  •  •  ,  Xn)  of 
the  same  degree  such  that  the  inequality 

(5.7)  Pl{Xl,  Xi,  “  ‘  ,  Xn)  ^  T[S{Xl,  X2,  -  Xn)]  ^  Piixi,  X2,  •  •  •  ,  Xn) 

is  satisfied  for  all  Xi,  X2,  -  ,  Xn.  A  quasipolynomial  is  said  to  be  admissible  if 

Pi  is  an  admissible  polynomial.  We  give  next  a  condition  which  assures  that 
fj{t)  is  the  characteristic  function  of  a  normal  distribution. 

Corollary  to  theorem  5.3.  Suppose  that  the  conditions  of  theorem  5.3  are 
satisfied  and  that  the  characteristic  function  fj{z)  has  no  zeros  in  the  entire  complex 
plane,  then  the  random  variable  Xj  is  normally  distributed. 

According  to  theorem  5.3,  the  function  f,{z)  is  an  entire  function  of  finite 
order  m  without  zeros.  We  apply  Hadamard’s  factorization  theorem  and  see 
th&t  fj{z)  =  exp  [Pm (2)]  where  Pm  is  a  polynomial  of  degree  m.  The  statement  of 
the  corollary  follows  then  from  Marcinkiewicz’  theorem. 

This  corollary  is  of  no  immediate  use  in  characterization  problems  since  one 
of  its  conditions  is  expressed  in  terms  of  the  characteristic  function  fj{t). 
It  is  therefore  desirable  to  find  a  condition  which  the  polynomial  statistic 
P(Xi,  •  •  •  ,  Xn)  of  theorem  5.3  must  satisfy  in  order  that  fj{z)  should  have  no 
(real  or  complex)  zeros.  Before  proceeding  further  we  must  introduce  a  special 
class  of  polynomials.  Let 

(5.8)  P(xi,  •  -  ,Xn)  =  L  Aj,...j„x{  ■  ■  •  xi^ 

Ji-i - H»  Sp 

be  a  polynomial  of  degree  p;  it  can  be  written  as  the  sum 

(5.9)  P(Xi,  •  •  •  ,  Xn)  =  Po(Xi,  •  •  •  ,  Xn)  +  Pl(Xi,  •  •  *  ,  Xn), 
where 

(5.10)  Po(xi,  •  •  •  ,  Xn)  =  L  x{‘  •  •  •  xi" 

- hy»=p 

is  a  homogeneous  polynomial  of  degree  p,  while  Pi(xi  •  •  •  ,  x„)  is  a  polynomial 
of  degree  less  than  p.  We  say  that  the  polynomial  P(xi,  •  •  •  ,  Xn)  is  nonsingular 
if  the  following  two  conditions  are  satisfied: 

(i)  Po(xi,  •  •  •  ,  Xn)  contains  the  pth  power  of  at  least  one  variable, 

(ii)  iro(p)  5*^  0  for  all  integers  p  >  0.  Here  7ro(p)  is  the  polynomial  formed  by 
replacing  each  positive  power  x|  by  =  p(p  —  1)  •  •  •  (p  —  j  +  1)  in 
Po(xi,  •  •  •  ,  Xn). 
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We  call  7ro(v)  the  adjoint  polynomial  to  P{xi,  •  •  •  ,  x„). 

Theorem  5.4.  Let  X\,  X2,  •  •  •  ,  X„  be  n  independently  and  identically  dis¬ 
tributed  random  variables  and  assume  that  the  characteristic  junction  j{z)  of  their 
common  distribution  is  an  entire  function.  Suppose  that  a  nonsingular  polynomial 
statistic 

(5.11)  P  =  P(Xi,  •  •  • ,  X„)  =  E  x^  •  -Xfr 

jl^ - 

of  degree  p  has  constant  regression  on  A  =  Xi  +  •  •  •  +  X„.  The  characteristic 
function  f{z)  has  then  no  zero  in  the  whole  complex  plane. 

We  emphasize  that  the  characteristic  function  f{t)  is  defined  for  all  complex 
arguments  z  =  t  -\-  iy,  where  t,  y  real,  by  writing  it  as  f{z).  Theorem  5.4  is  due 
to  Yu.  V.  Linnik  [29] ;  we  give  here  a  somewhat  simplified  proof. 

It  easily  follows  from  lemma  4.1  that  the  relation 

(5.12)  X(Pe*^A)  =  EiP) 
holds  for  all  complex  z.  We  write 

(5.13)  /<'■)  =  /(J)(z)  =  =  ii  E(Xie<’^) 


and  note  that  f^^''(z)  =  f(z).  We  see  from  (5.12)  that 


(5.14)  E  ,--(y.+  ---+A)  .  ./O")  =  E{P)[m]\ 

AH - 


We  give  an  indirect  proof  of  the  theorem  and  assume  therefore  that  the  func- 
tion/(2)  has  zeros.  Let  the  point  z  =  s’o  be  one  of  the  zeros  oif{z)  which  are  nearest 
to  the  origin  and  denote  the  order  of  z^  by  v,  where  p  is  a  positive  integer. 

Since  /(z)  does  not  vanish  in  the  circle  lz|  <  Izoj  we  may  divide  (5.14)  by 
[/(z)]".  We  write 

(5.15)  7?o(z)  =  E 

jl-] - \-jn  =  p  J 


and 

(5.16) 


/i'i(z)  =  E 

ji-\ - hin  <P 


fUO  .  .  .  f(u) 


and  see  that 


(5.17)  Poiz)  +  Rl{z)  —  ;zl  <  Izoi, 

where  C  =  i’E{P).  Let 

(5.18)  <#>  =  <t>iz)  =  log/(2); 
it  is  then  easily  verified  that 

(5.19)  -y  =  +  eM',  (/)",  •  •  •  ,  <#>('■-')),  i  =  1,  2,  •  •  •  , 

where  dj  is  a  polynomial  in  <(>',  <t>",  •  •  •  ,  We  also  write  =  1  and  do  =  0. 

We  substitute  (5.19)  into  (5.17)  and  get  for  [z]  <  Izo] 


(5.20) 


iS:«(z)  +  Sdz)  =  C, 
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where 


(5.21) 

and 


«o(2)  =  L  Aj...  +  9^,] .  .  . 

AH - 1-71.  =  p 


(5.22)  S.(2)  =  £ 

AH — i-yi.<p 
According  to  our  assumptions 

(5.23)  /(2)  =  (z  -  2<,)V(«), 

where  ^(z)  is  an  entire  function  such  that  g(zii)  0.  It  is  easy  to  verify  that 


(5.24)  0'(z)  =  +  5,(2), 

Z  —  Zo 

and  in  general 

(5.25)  4,«i(2)  =  +  5,(2),  i  =  1,  2,  •  •  •  . 


Here  the  functions  h,{z)  are  regular  at  2  =  Zq.  We  substitute  (5.25)  into  (5.20) 
and  see  that 


(5.26) 


_ l£ _  4.  _ l£z:l _ 4.  ...  4  -J}  . 

(2  —  Zo)p  (z  —  Zo)^^  z  —  Zo 


+  H(Z)  =  C, 


where  ^(2)  is  regular  at  2  =  20. 

We  note  that  relation  (5.26)  leads  to  a  contradiction  if  at  least  one  of  the 
coefficients  71,  72,  •  •  •  ,  7p  is  different  from  zero.  We  complete  the  proof  of  the 
theorem  by  showing  that  jp  9^  0. 

We  remark  that  jp  depends  only  on  v  and  on  the  coefficients  of  the  homo¬ 
geneous  polynomial  Poixi,  X2,  •  •  •  ,  Xn).  We  see  that  jp  is  the  coefficient  of 
(2  —  2o)~P  in  the  expression  which  we  obtain  by  substituting  (5.25)  into  (5.21). 
We  get  the  same  value  for  the  coefficient  of  (2  —  zo)~p  if  we  substitute  (5.23)  into 
(5.15).  The  coefficient  7p  can  also  be  obtained  by  substituting  ^(2)  =  Ci{z  —  20)” 
instead  of  f{z)  into  (5.15).  Here  (7i  0  is  a  constant.  We  see  then  that 


(5.27) 

Therefore 

T 

I 

II 

(5.28) 

7p  =  L  •  • 

AH — i-y»=p 

Thus  7p  =  Tro{v)  7^  0  for  all  positive  integer  values  of  v  and  theorem  5.4  is 
proved.  We  can  use  theorems  5.3  and  5.4  and  the  corollary  to  theorem  5.4  to 
get  the  following  characterization  of  the  normal  distribution. 

Theorem  5.5.  Let  Xi,  X2,  •  •  •  ,  Xn  he  n  independently  and  identically  dis¬ 
tributed  random  variables.  Let  P  =  P{Xi,  X2,  •  •  •  ,  X„)  he  an  admissible,  non¬ 
singular  polynomial  statistic  and  A  =  Ai  -f  X2  +  •  •  •  +  Xn.  If  P  and  A  are 
independent  then  the  common  distribution  of  the  random  variables  Xj  is  normal. 

We  consider  apphcations  of  theorem  5.5. 

Theorem  5.6.  Let  Xi,  X2,  •  •  •  ,  Xn  he  n  independently  and  identically  distri- 
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huted  random  variables  with  common  distribution  function  F(x).  Let  A  =  Xi  + 
X2  +  •  •  •  +  and  P  =  P{Xi,  X2,  •  •  •  ,  Xn)  be  an  admissible,  homogeneous  poly¬ 
nomial  statistic  of  degree  p  satisfying 

(i)  P  and  A  are  independently  distributed, 

(ii)  E{P)  =  Kp  where  Kp  is  the  pth  cumulant  of  F(_x).  Then  F{x)  is  normal. 

Note  that  the  existence  of  the  cumulants  of  F{x)  is  insured  by  the  fact  that 

F{x)  has,  according  to  theorem  5.4,  an  entire  characteristic  function.  Let  f{t)  be 
the  characteristic  function  of  F(x).  Then  <l>(t)  —  log  f(t)  exists  in  a  certain  neigh¬ 
borhood  of  the  origin.  It  follows  from  the  assumptions  of  our  theorem  that 
^(p)(^)  =  dp<l>(t)/dt^  exists.  It  is  easily  seen  by  induction  that  <l>^^^(t)  is  a  poly¬ 
nomial  in  /'//,  /"//,  •  •  •  ,  /^^V/  which  has  the  form 


The  summation  is  here  extended  over  all  nonnegative  integers  Si,  sz,  ■  •  •  ,  Sp 
which  satisfy  the  relation 

(5.30)  Si  -}-  2s2  +  •  •  •  pSp  =  p. 

We  put  <  =  0  in  (5.29)  and  obtain 

(5.31)  Up  =  H  haist-  •  spOi^ioi2-  •  -ap  =  G{ai,  0:2,  •  •  •  ,  (Xp), 

where  aj  is  the  jth  moment  of  F(x)  and  where  the  summation  is  extended  over 
all  subscripts  satisfying  (5.30).  Let 

(5.32)  P=  L 

ji-l - h/«  =  p 

it  is  always  possible  to  choose  coefficients  so  that 

(5.33)  P  = 

Here,  and  in  the  following,  the  summation  ^  is  extended  over  all  Si,  •  •  •  ,  Sp 
satisfying  (5.30)  while  the  summation  runs  over  all  permutations  of 
Uh  •  •  •  >  fn)  such  that  Sr  of  these  exponents  equal  r  for  r  =  1,2,  •••,?>  and  the 
remaining  n  —  Si  —  •  •  •  —  Sp  exponents  are  zero. 

It  follows  from  (5.33)  that 

(5.34)  P(P)  =  L  K-sp  •  -4”  E*  5*... A. 

Hence  we  see  from  condition  (i)  and  (5.31)  that 

(5.35)  E*  •  •  •»,;  n  •  •  -in  =  1 

for  all  p-tuples  (si,  •  •  •  ,  Sp)  for  which  X*...,,  7^  0.  We  form  now  the  polynomial 
To(v)  corresponding  to  P.  We  get  from  (5.33) 

(5.36)  iro(v)  =  E  X,i....p 
and  thus  in  view  of  (5.35)  and  (5.31) 

(5.37)  7ro(v)  =(7(p<»,  •••  ,p(^>). 

Since  Kp  is  the  coefficient  of  t^/p\  in  the  expansion  of  log  [Ey  ay^ViHj  we  note 
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that  ro(i')  is  the  coefficient  of  t^/pl  in  the  expansion  of  log  = 

iog  [(1  +  O']-  We  obtain  therefore 

(5.38)  iro(v)  =  (  — —  1)! 

Since  xo(j')  9^  0  for  all  positive  integers  v,  we  see  that  P  is  a  nonsingular  poly¬ 
nomial.  Theorem  5.6  follows  then  immediately  from  theorem  5.2. 

The  conditions  of  theorem  5.6  are  satisfied  if  P  is  the  A:-statistic  of  order  p. 
This  special  case  was  already  proved  earlier  by  several  authors  (see  [38],  p.  201). 
However,  it  was  then  assumed  that  the  pth  moment  of  F{x)  exists.  The  more 
general  methods  used  here  make  this  assumption  unnecessary. 

A  second  class  of  polynomial  statistics  for  which  it  is  possible  to  compute  the 
adjoint  polynomial  are  the  central  moments. 

Lemma  5.1.  Let  p  be  a  positive  integer  and  write 

(5.39)  P  =  PiXi,  •  •  •  ,  .Tn)  =  Z  (xk  —  x)^,  X  =  -  X]  Xic. 

k=i  nk=i 


The  adjoint  polynomial  of  P  is  then 


(5.40) 


TTpiv) 


(-1)^7?!  f 
nP-^  r=o 


(-mn  -  1) 


For  the  proof  of  this  lemma  we  refer  to  [22].  We  can  use  the  explicit  form  of 
Tp{v)  to  derive  a  condition  which  assures  that  the  sample  moment  of  order  p  is 
a  nonsingular  polynomial  statistic. 

Lemma  5.2.  Let  pbea  positive  integer  and  write 


(5.41) 


P  =  P{Xi,  •  "  ,Xn) 


=  L  (xk 

k=i 


x)p, 


1  ” 

x  =  -Zx,. 
nk=i 


If  {p  —  1)!  is  not  divisible  by  n  —  I,  the  adjoint  polynomial  Tp{v)  of  P  has  no 
nonzero  integer  roots. 

Suppose  that  for  some  integer  v,  with  v  9^  0,  we  have  xp(j')  =  0.  Then  we  see 
from  lemma  5.1  that 

- 1"' 

Thus  multiplying  by  pi  and  cancelling  the  common  factor  (n  —  I)*/  we  find  that 
(5.43)  (nv  —  V  —  l)(nv  —  p  —  2)-  ■  -{nv  —  v  —  p  -]r  L)  =  0{mod  {n  —  1)} 


so  that  , 

(5.44)  (p  —  1)!  =  0{mod  (n  —  1)}. 

We  consider  in  the  following  a  sample  X\,  Xt,  •  •  •  ,  X„  from  a  certain  population 
and  write 

(5.45)  J  =  -{X,  + X,+  ■■■  +  X,),  m,  =  i  t  (X;  -  X)- 

n  n 
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for  the  sample  mean  and  the  sample  central  moment  of  order  p,  respectively. 
We  prove 

Theorem  5.7.  Let  X\,  X<i,  •  •  •  ,  X„  he  a  sample  of  size  nfrom  a  certain  popula¬ 
tion.  Let  pbe  a  positive  integer  such  that  (p  —  1) !  is  not  divisible  by  (n  —  1).  The 
population  is  normal  if,  and  only  if,  the  sample  central  moment  mp  of  order  p  is 
independently  distributed  of  the  sample  mean  X. 

Remark.  The  condition  that  (p  —  1) !  is  not  divisible  by  (n  —  1)  is  satisfied 
if  w  >(p  -  1)!  +  1. 

The  necessity  of  the  condition  of  theorem  5.7  follows  from  the  well  known  fact 
that  in  a  normal  population  any  translation-invariant  statistic  is  independent 
of  the  sample  mean.  The  sufficiency  of  the  condition  follows  from  theorem  5.5  and 
lemma  5.2. 

We  conclude  this  section  by  mentioning  two  other  recent  developments. 

Yu.  V.  Linnik  [28]  studied  the  possibility  of  determining  the  family  of  dis¬ 
tribution  functions  to  which  the  population  distribution  function  belongs  from 
the  distribution  function  of  a  statistic.  The  second  problem  is  of  a  different 
nature  since  it  concerns  the  characterization  of  distributions  and  not  of  popula¬ 
tions.  Let  X  and  Y  be  two  independently  and  identically  distributed  random 
variables  and  assume  that  their  common  distribution  is  normal  with  zero  mean. 
It  is  then  known  that  the  quotient  Z  =  X/Y  is  distributed  according  to  a  Cauchy 
law.  A  number  of  authors  investigated  whether  the  normal  distribution  could  be 
characterized  by  this  property.  Counterexamples  were  constructed  by  R.  G.  Laha 
[15],  G.  P.  Steck  [56],  and  J.  G.  Mauldon  [42].  R.  G.  Laha  [16],  [17]  undertook 
a  systematic  study  of  the  family  C  of  distributions  F{x)  which  have  the  following 
property :  If  X  and  Y  are  two  independently  and  identically  distributed  random 
variables  with  common  distribution  F{x),  then  the  quotient  X/Y  has  a  Cauchy 
distribution.  He  obtained  a  number  of  properties  of  F{x)  and  could  also  charac¬ 
terize  the  normal  distribution,  assuming  F(x)  G  C  and  some  additional  restric¬ 
tions  which  include  the  existence  of  all  moments  of  F(x) . 


Part  III.  The  arithmetic  of  distribution  functions  and  related  topics 

The  arithmetic  of  distribution  functions  deals  with  the  decomposition  of 
characteristic  functions  into  factors  which  are  characteristic  functions  of  non¬ 
degenerate  distributions.  Most  of  the  recent  developments  in  this  area  were 
stimulated  by  some  results  obtained  by  H.  Cramer  [4],  P.  L6vy  [24]  and 
D.  A.  Raikov  [47]  more  than  twenty  years  ago.  We  state  here  two  of  the  classical 
results. 

Cramer’s  theorem.  Letf{t)  =  exp  {iyt  —  (tH'^12]  he  the  characteristic  function 
of  the  normal  distribution  and  suppose  that  f{t)  =  fi{t)f2{t)  where  fi(t)  and  f2(t) 
are  characteristic  functions.  Then  /i(0  and  f2{t)  are  necessarily  characteristic  func¬ 
tions  of  normal  distributions. 

This  theorem  was  first  conjectured  by  P.  L6vy  and  later  proved  by  H.  Cram6r 
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and  it  is  therefore  often  called  the  theorem  of  L4vy-Cram4r.  A  similar  property 
of  the  Poisson  distribution  was  established  by  D.  A.  Raikov. 

Raikov’s  THEOREM.  Letf{t)  =  exp  [X(e’‘  —  1)]  he  the  characteristic  function  of 
the  Poisson  distribution  and  suppose  that  f(t)  =  fi(t)f2{t)  where  fi(t)  and  f2{t)  are 
characteristic  functions.  Then  fi(t)  andf2{t)  are  necessarily  characteristic  functions 
of  Poisson  distributions. 

We  see  immediately  from  the  theorems  of  Cramer  and  Raikov  that  the  normal 
distribution,  as  well  as  the  Poisson  distribution,  belongs  to  the  family  of  infinitely 
divisible  distributions  that  have  no  indecomposable  factors.  We  shall  discuss 
this  family  in  section  6,  and  we  shall  treat  other  problems  which  are  closely 
connected  with  the  investigations  of  Cram4r,  L4vy,  and  Raikov  in  sections  7 
and  8. 

6.  Infinitely  divisible  characteristic  functions  that  have  no 

indecomposable  factors 

D.  A.  Raikov  raised  (in  [47])  several  questions  concerning  the  structure  of 
infinitely  divisible  characteristic  functions  that  have  no  indecomposable  factors. 
These  problems  arc  very  difficult  and  the  first  advances  in  this  area  were  made 
by  Yu.  V.  Linnik  approximately  twenty  years  after  the  publication  of  Raikov’s 
paper.  Yu.  V.  Linnik’s  work  on  the  factorization  of  infinitely  divisible  laws  will 
be  discussed  in  this  section. 

D.  A.  Raikov  and  P.  L4vy  established  the  following  remarkable  fact:  The 
convolution  of  two  Poisson  type  distributions  has  no  indecomposable  factors; 
however,  convolutions  of  three  Poisson  type  distributions  may  have  indecom¬ 
posable  factors.  This  result  made  it  desirable  to  investigate  the  factorization  of 
the  convolution  of  a  normal  and  a  Poisson  distribution.  Yu.  V.  Linnik  [30],  [31] 
obtained 

Theorem  6.1.  Let 

(6.1)  f(t)  =  exp  *fx(e‘‘  —  1)  +  int  —  1’  M  real,  a'  ^0,  X  ^  0, 


be  the  characteristic  function  of  the  convolution  of  a  normal  and  of  a  Poisson  distri¬ 
bution.  Suppose  thatfit)  admits  the  decomposition  f{l)  =  /i (0/2(0-  Then 


(6.2) 

fj{t)  =  exp  |xXe*‘ 

-  1)  +  - 1  ,.5(4- 

j  =  1,  2, 

where 

(6.3) 

X  =  Xi  -b  X2, 

o'  =  Ol  +  O2,  Xy  ^  0, 

0 

All 

b 

Theorem  6.1  contains  Cramer’s  theorem  and  Raikov’s  theorem  as  special 
cases.  However,  the  method  of  proof  is  entirely  different  and  requires  more 
powerful  anal3rtical  tools.  This  is  explained  by  the  fact  that  Cramer’s  theorem 
deals  with  an  entire  characteristic  function  of  finite  order  while  Raikov’s  theorem 
is  concerned  with  a  periodic  characteristic  function  with  real  period  (thus  nec¬ 
essarily  the  characteristic  function  of  a  lattice  distribution).  Under  the  assump- 
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tions  of  theorem  6.1  both  these  advantages  are  lost,  and  the  proof  becomes  much 
more  complicated  and  requires  some  sharp  estimates  used  in  Vinogradov’s  study 
of  trigonometric  sums. 

The  method  developed  for  the  proof  of  theorem  6.1  also  permitted  the  attack 
of  the  more  general  factorization  problems  of  infinitely  divisible  laws.  A  series 
of  papers  by  Yu.  V.  Linnik  [32],  [33],  [34],  [35],  [36]  contains  his  investigations 
on  this  subject;  the  principal  problem  is  the  determination  of  the  structure  of  the 
infinitely  divisible  laws  that  have  no  indecomposable  factors.  In  the  following, 
we  call  this  class  of  characteristic  functions  the  class  /q.  In  order  to  discuss 
Linnik’s  investigations  of  the  class  /o  we  must  introduce  certain  terms  and 
notation. 

It  is  well  known  that  every  infinitely  divisible  characteristic  function  f(t)  can 
be  written  in  the  canonical  form 

(6.4)  log  f(t)  =  ita  —  +  j  -  1  - 

where  the  constants  a  and  are  real,  ^  0,  and  where  M{u)  and  N{u)  satisfy 

(i)  M{u)  and  N(u)  are  nondecreasing  in  the  intervals  (  — «,  0)  and  (0,  +») 
respectively, 

(ii)  the  integrals  dM{u)  and  j*  dN{u)  are  finite  for  every  e  >  0, 

(iii)  M{  —  oo)  =  =  0. 

It  is  convenient  to  call  M{u)  the  negative  and  N(u)  the  positive  Poisson 
spectrum  of  f(t).  We  say  that  an  infinitely  divisible  characteristic  function  has 
bounded  negative  Poisson  spectrum  if  there  exists  a  positive  number  A  such 
_ ^ 

dM{u)  =  0.  Similarly  we  define  infinitely  divisible  distributions  with 

bounded  positive  Poisson  spectrum.  We  say  that  the  Poisson  spectrum  of  an 
infinitely  divisible  characteristic  fimction  f(t)  is  bounded,  if  its  positive  as  well 
as  its  negative  Poisson  spectrum  is  bounded.  An  infinitely  divisible  character¬ 
istic  function  is  said  to  have  a  finite  spectrum  if 

1  w  n 

(6.5)  log/(<)  =  ait  -  5  aV  +  Z  -  D  +  Z  -  1). 

^  y-1  y-i 

Here  m  and  n  are  nonnegative  integers,  Xy  >  0,  X-y  >  0,  My  >  0,  vj  >  0.  If 
either  m  orn  is  equal  to  zero  then  the  corresponding  sum  is  omitted.  An  infinitely 
divisible  characteristic  function /(/)  is  said  to  have  a  denumerable  Poisson  spec¬ 
trum  if 

(6.6)  log  fit) 

=  aU  -  I  a’e  +  ,z  \  («“"  -  1  -  f^)  +  t  X-,  (e-"-  -  1  +  Y^). 
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where  \j  >  0,  >  0  and  where  the  series 


(6.7) 


y  =  l  1  + 


and 


i=i  1  + 


are  convergent  and  where 

(6.8)  E  +  E 

/i)  <e  3>i  <t 

tends  to  zero  as  e  approaches  zero.  The  numbers  Mi  and  vj  are  called  the  Poisson 
frequencies,  the  Xy  and  X_,  the  energy  parameters  of  f{t).  A  characteristic  function 
with  finite  or  denumerable  Poisson  spectrum  is  said  to  have  a  rational  spectrum 
if  nj/nk  =  Tjk  and  Vj/pk  =  sjk  are  rational  numbers.  Let/(0  be  a  characteristic 
function  with  a  bounded  Poisson  spectrum  which  is  contained  in  the  interval 
We  say  that  the  spectrum  of  f(t)  is  rational  to  the  right  of  the  point 
a,  where  0  <  a  <  if 

(6.9)  log  m  =iia-\  +  j  ■“  -  1  -  j^,)  dM(u) 

+ /;  (^‘'“  - 1  -  rli?)  h  (^)  - 1] 

where  m  >  0,  X,  >  0,  where  y  =  1,  •  •  •  ,  m,  while  0  <  ai  <  •  •  •  <  a„  are  integers 
such  that  ain/q  >  a.  In  a  similar  manner  we  can  define  a  Poisson  spectrum 
which  is  rational  to  the  left  of  the  point  where  —  A  <  /?  <  0. 

We  can  now  state  some  of  Linnik’s  results. 

Theorem  6.2.  In  order  that  an  infinitely  divisible  characteristic  function  with 
normal  component,  with  >  0,  should  have  no  indecomposable  factor  it  is  necessary 
that  its  Poisson  spectrum  be  finite  or  denumerable.  Moreover,  the  Poisson  frequencies 
of  the  positive  spectrum  must  have  the  form 


(6.10) 


while  the  Poisson  frequencies  of  the  negative  spectrum  must  be 


(6.11) 


•  •  •  ,  g-zg-ip,  g-iv, 


V  V 

V,  —J  - > 

g\  g-iQt 


where  the  •  •  • ,  /c_2,  /c-i,  ki,  /c2,  •  •  •  and  the  •  •  • ,  gr_2,  g-i,  gi,  g^,  •  •  •  are  arbitrary 
integers  {not  necessarily  all  different)  which  are  greater  than  one.  If  the  Poisson 
spectrum  of  f{t)  is  bounded  then  this  condition  is  also  sufficient  for  the  absence  of 
indecomposable  factors. 

This  is  probably  one  of  the  most  important  advances  in  the  arithmetic  of 
distribution  function  since  the  publication  of  the  fundamental  studies  of  P.  L4vy, 
A.  I.  Khinchin,  and  D.  A.  Raikov.  The  necessity  of  the  condition  of  theorem  6.2 
is  established  in  [34].  The  proof  is  based  on  several  lemmas  which  are  of  inde¬ 
pendent  interest.  As  an  example  we  mention  the  following  theorem  (this  is 
Linnik’s  basic  lemma  I). 
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Theorem  6.3.  Let  7,  Xi  and  X2  he  three  positive  numbers  and  let  a  =  p/q  he  a 
rational  number  where  p  and  q  are  two  integers  such  that  I  <  p  <  q  and  (p,  5)  =  1 . 
Then 


(6.12)  f{t)  =  exp  [— 7^^  +  Xi(e^‘  —  1)  +  —  1)  —  —  1)] 

is  a  characteristic  function  provided  that  v  >  0  is  sufficiently  small. 

Remark.  This  theorem  can  be  used  to  construct  convolutions  of  a  normal, 
a  Poisson,  and  a  Poisson  type  distribution  that  have  indecomposable  factors. 
This  is  similar  to  the  result  of  D.  A.  Raikov  and  P.  L6vy  concerning  the  factors 
of  a  convolution  of  three  Poisson  type  distributions  (see  the  paragraph  before 
theorem  6.1). 

The  sufficiency  of  the  condition  of  theorem  6.2  for  characteristic  functions 
with  bounded  Poisson  spectrum  is  established  in  [35].  In  the  last  mentioned 
paper  we  find  also  the  proofs  of  several  theorems,  stated  already  in  [34].  We 
quote  here  only  one  of  these  results. 

Theorem  6.4.  Let  f(t)  he  a  characteristic  function  with  hounded  Poisson  spec¬ 
trum,  then  all  its  factors  have  the  form 


(6.13) 


g(t)  =  exp  {Pziit)  +  t*  ^(li)  du}, 


where  Pz{it)  is  a  polynomial  of  degree  not  exceeding  three  and  where  ^{u)  is  quad- 
ratically  integrahle  (in  the  Lehesgue  sense)  over  the  interval  [—A,B]  which  contains 
the  spectrum  of  f(t). 

Linnik  [35]  studied  infinitely  divisible  characteristic  functions  that  have  a 
normal  component  and  a  bounded  Poisson  spectrum  which  is  rational  to  the 
right  (or  left)  of  a  point  of  [—A,  B].  He  determined  the  possible  factors  of  such 
characteristic  functions. 

The  investigation  of  infinitely  divisible  characteristic  functions  without  a  nor¬ 
mal  component  or  of  infinitely  divisible  characteristic  functions  with  unbounded 
Poisson  spectrum  is  very  difficult.  We  mention  here  one  of  the  few  results  (see 
[36])  which  are  known  at  present. 

Theorem  6.5.  Suppose  that  an  infinitely  divisible  characteristic  function  f(t) 
has  the  form  (6.6)  and  that  its  Poisson  frequencies  m  and  Vj  satisfy  the  conditions 
(6.10)  and  (6.11).  Assume  further  that  there  exist  positive  constants  n,  v,  c  and  a  such 
that 


(6.14) 


log  log  ^  >  C 
log  log  ^  >  c  v] 


for  all  m  >  n  and  pj  >  v,  then  f(t)  has  only  infinitely  divisible  components. 

The  condition  of  the  theorem  means  that  the  energy  parameters  corresponding 
to  high  frequencies  decrease  rapidly. 

The  theorems  discussed  in  this  section  indicate  that  the  class  /o  forms  a  rather 
small  subset  of  the  class  of  all  infinitely  divisible  characteristic  functions.  The 
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problem  of  finding  the  indecomposable  factors  of  infinitely  divisible  characteris¬ 
tic  functions  not  belonging  to  the  class  h  has  not  yet  been  attacked.  For  instance, 
it  is  known  (as  a  consequence  of  a  more  general  result  of  H.  Cramer  [5])  that 
the  Gamma  distribution  has  indecomposable  factors;  however,  these  factors  have 
not  yet  been  determined. 


7.  Factor-closed  and  strongly  factor-closed  families 

The  theorems  of  Cramer  and  Raikov  which  we  stated  in  the  introductory 
paragraph  of  part  III  can  be  formulated  in  a  somewhat  different  manner. 

We  say  (see  H.  Teicher  [57])  that  a  family  f  of  characteristic  functions  is 
factor-closed  if  the  relations 

(i)  /  G  f 

(ii)  /  =  /1/2,  where  /i,  /a  are  characteristic  functions,  imply  that  fi  G  f,  /2  G  f • 

The  theorems  of  Cram4r  and  Raikov  admit  then  the  following  two,  different 

formulations:  (a)  The  characteristic  functions  of  normal  (respectively  Poisson) 
distributions  belong  to  the  class  /o.  (b)  The  characteristic  functions  of  normal 
(respectively  Poisson)  distributions  form  a  factor-closed  family. 

The  first  formulation  leads  to  the  studies  discussed  in  section  6;  in  the  following 
we  treat  problems  which  are  related  to  the  second  formulation. 

A.  A.  Zinger  and  Yu.  V.  Linnik  [65]  stimulated  these  investigations  by  proving 
the  following,  interesting  extension  of  Cramer’s  theorem. 

Theorem  7.1.  Letfi(t),  •  •  •  ,  fn(t)  be  arbitrary  characteristic  junctions  and 
suppose  that  the  relation 

(7.1)  n  [/i(0]“’'  =  exp  [int  -  ^  <tH^] 

holds  for  some  positive  real  numbers  ai,  a2,  •  •  •  ,  ««  in  some  neighborhood  of  the 
origin.  Then  the  characteristic  functions  f  j{t) ,  for  j  =  1,  2,  •  •  •  ,  n,  belong  to  normal 
distributions. 

In  order  to  give  a  concise  formulation  of  certain  results  we  introduce  the  fol¬ 
lowing  terminology. 

A  family  f  of  characteristic  functions  is  said  to  be  strongly  factor-closed  if 
the  relations 

(i)  /Gf 

(ii)  IT-  i/P  =  f,  where /i,  •  •  •  ,  /n  are  characteristic  functions,  valid  for  some 
positive  ai,  •  •  •  ,  a„  imply  that/y  G  f,  forj  =  1,  2,  •  •  •  ,n. 

We  can  then  state  the  result  of  theorem  7.1  by  saying  that  the  family  of 
normal  distributions  is  strongly  factor-closed. 

Every  strongly  factor-closed  family  is  also  factor-closed;  at  present  one  has  no 
example  of  a  factor-closed  family  of  characteristic  functions  which  is  not  strongly 
factor-closed.  The  properties  of  being  factor-closed  or  strongly  factor-closed  have 
a  different  character:  The  first  can  be  regarded  as  a  probabilistic  property  of  a 
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family  of  characteristic  functions  and  can  also  be  expressed  in  terms  of  random 
variables;  the  second  is  essentially  of  an  analytical  nature.  It  is  therefore  not 
surprising  that  the  proofs  establishing  these  two  properties  use  different  tech¬ 
niques. 

In  conclusion,  we  give  a  table  of  references  to  papers  which  contain  the  proofs 
that  certain  families  of  characteristic  functions  are  factor-closed,  respectively, 
strongly  factor-closed.  In  this  table  we  denote  the  family  of  analytic  character¬ 
istic  functions  by  f^,  the  family  of  entire  characteristic  functions  and  write 
9K2OT  for  the  family  of  characteristic  functions  which  have  derivatives  up  to  the 
order  2m,  where  m  is  a  positive  integer.  Actually,  H.  Teicher  [57]  discusses  a 
somewhat  more  general  family  than  the  binomial. 


TABLE  II 


Family 

Factor-closed 

Strongly  factor-closed 

Normal 

Cram4r  [4] 

Zinger-Linnik  [65] 

Poisson 

Baikov  [47] 

Dugu6  [7],  [8] 

Binomial 

Teicher  [57] 

Teicher  [58] 

f.1. 

Baikov  [47], 

L6vy  [23],  [25] 

Laha  [14] 

Devinatz  [6] 

proof  similar  to  the  one  in 
[65]  makes  use  also  of  [6] 

We  mention  also  that  the  family  of  the  characteristic  functions  of  discrete 
distributions  is  factor-closed  and  that  the  characteristic  functions  of  lattice  dis¬ 
tributions  form  a  strongly  factor-closed  family.  The  proofs  of  these  two  state¬ 
ments  are  almost  trivial.  Yu.  V.  Linnik  [37]  showed  also  that  a  rather  wide 
subset  of  the  class  /o  is  strongly  factor-closed. 

Yu.  V.  Linnik  [27]  questioned  whether  it  is  possible  to  extend  theorem  7.1  to 
infinite  products.  This  problem  was  solved  recently  by  L.  V.  Mamay  [41]  who 
obtained 

Theorem  7.2.  Let  for  j  =  1,2,  •  he  an  arbitrary  sequence  of  char¬ 

acteristic  functions  and  let  {affjfor  f  =  1,  2,  •  •  •  ,hea  sequence  of  positive  numbers 
such  that  aj  >  €0  >  0.  Suppose  that  g{z)  is  a  function  of  the  complex  variable 
z  =  t-\-  iy  which  is  regular  in  the  strip 

(7.2)  llm(2)l  <  X 

and  which  has  no  zeros.  Assume  further  that  there  exists  a  A  >  0  such  that 

(7.3)  n  =  gil) 

J-1 

is  valid  for  <  A.  Then  the  fj{t)  are  analytic  characteristic  functions  which  are 
regular  at  least  in  the  strip  (7.2)  and  the  relation  (7.3)  hclds  also  for  complex  values 
satisfying  (7.2). 
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8.  Stability  theorems 

We  say  that  a  family  of  distribution  functions  is  factor-closed  if  the  cor¬ 
responding  characteristic  functions  form  a  factor-closed  family.  The  existence  of 
factor-closed  families  leads  to  the  following  question.  Suppose  that  a  distribution 
function  F{x)  is,  in  some  sense,  close  to  a  distribution  of  a  factor-closed  family  f 
and  assume  that  F  is  the  convolution  of  tw’o  distributions  Fi  and  F^,  that  is, 
F  =  Fi*Fi.  Is  it  then  possible  to  assert  that  the  components  Fi  and  F2  are  neces¬ 
sarily  close  to  some  distribution  of  the  factor-closed  family  f?  We  call  theorems 
which  make  this  assertion  stability  theorems. 

In  the  following  we  denote  the  normal  distribution  function  by 

(8.1)  ^{x)  = j 

and  the  Poisson  distribution  function  by 

(8.2)  F(x-,\) 

Here  e{x)  =  0  for  a;  <  0  but  €(x)  =  1  for  rc  ^  0  is  the  degenerate  distribution 
function,  while  X  is  a  positive  constant. 

The  first  stability  theorem  was  derived  by  N.  A.  Sapogov  [48]  who  obtained 
the  following  result. 

Theorem  8.1.  Let  X  =  Xi  X2  be  the  sum  of  two  independent  random  var¬ 
iables  and  assume  that  the  distribution  function  F{x)  of  X  satisfies  the  condition 

(8.3)  sup  lP(a;)  —  4>(a;)l  <  €, 

— «  <*  <  * 

where  e  <  1  is  a  given  positive  number.  Let  Fi{x)  be  the  distribution  function  of  Xi 
and  write 

xdFi{x),  ^  I  ^^dFi(x)  —  ai  >  0,  iV  =  yog-J  . 

Then 

(8.5)  sup  Fi{x)—^(- - (log ; 

here  C  is  a  constant  which  is  independent  of  e,  <ri,  ai.  A  similar  inequality  holds  for 
the  distribution  function  F2{x)  of  X2. 

Recently  N.  A.  Sapogov  [49],  [50]  improved  this  estimate  and  extended  it  also 
to  random  vectors.  We  state  also  a  stability  theorem  for  the  family  of  Poisson 
type  distributions  which  is  due  to  O.  V,  Shalayevskiy  [51]. 

Theorem  8.2.  Let  X  =  Xi  A-  X2  be  the  sum  of  two  independent  random  var¬ 
iables  and  assume  that  the  distribution  function  F{x)  of  X  satisfies  the  condition 

(8.6)  sup  |F(x)  —  F(x;X)|  <  e, 

—  00  <x  <  00 

where  e  <  1  and  X  are  positive  numbers.  Let  F,{x)  be  the  distribution  function  of 
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^31  for  j  =  1,2,  and  let  a  he  the  upper  hound  of  those  values  y  for  which  P{Xi  <  y} 

/—  rN+i  rN+i 

^  V  6.  We  write  L  ^  dFi{x  +  o)  and  ^2  =  L  ^  dF^ix  +  a),  where 
1  /e  =  N^.  Then  for  sufficiently  small  e 

sup  |Fi(a:)  -  F(a;  -  a;  Xi)]  <  (\  +  log  -1  , 

—«<*<«>  \  ^/L 

sup^  li^’aCa:)  -  F(x  +  a;  X2)  |  <  ^X  +  “]  > 

where  <a  <  1/2  is  a  constant. 

A  general  stability  theorem  for  a  wide  class  of  infinitely  divisible  distributions 
was  given  by  Yu.  V.  Linnik  [36]. 
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